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Abstract 

Applying a Lindemann-like criterion obtained previously by Kierfeld, Natter- 
mann and Hwa [Phys. Rev. B 55, 626 (1997)], we estimate the magnetic field 
and temperature for a high-T c superconductor, at which a topologically or- 
dered vortex glass phase becomes unstable with respect to a disorder-induced 
formation of dislocations. The employed criterion is shown to be equivalent 
to a conventional phenomenological Lindemann criterion including the values 
for the numerical factors, i.e., for the Lindemann-number. The positional 
correlation length of the topologically ordered vortex glass is calculated. 

I. INTRODUCTION 

The influence of disorder on the Abrikosov vortex lattice in the mixed phase of high- 
temperature superconductors, such as YBa 2 Cu 3 07_5 (YBCO), is an issue of immediate 
technological interest because pinning of the flux lines by disorder opens the possibility of 
regaining a dissipation-free current flow in the mixed phase. The flux line (FL) array in a 
high-temperature superconductor (HTSC) is extremely susceptible to thermal and disorder- 
induced fluctuations due to the interplay of several parameters, namely the high transition 
temperature T c , large magnetic penetration depths A and short coherence lengths £, and 
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a strong anisotropy of the material. This leads to the existence of a variety of fluctuation 
dominated phases of the FL array and very rich phase diagrams for the HTSC materials [|TJ. 
We want to consider here the pinning of FLs by point defects such as the oxygen vacancies, 
which is usually referred to as point disorder. It is well-known that the FL lattice is unstable 
to point disorder [Q. It has been conjectured that due to a collective pinning by the point 
disorder, the FL array may form a vortex glass phase with zero linear resistivity || [| . 
Although the existence of a vortex glass (VG) phase has been verified experimentally fl7HTT 



its large scale properties characterizing the nature of the VG phase are still under debate [|TJ . 
A possible scenario for a description of the low-temperature properties of the FL array 
subject to point disorder is the existence of a topologically ordered, i.e., dislocation-free 
VG phase, the so-called Bragg glass phase |12| as a thermodynamically stable phase. In this 



glassy phase, a quasi long range positional order of the FL array is maintained in spite of the 
pinning P,[12|. This entails the existence of algebraically decaying Bragg peaks in diffraction 
experiments on this phase (bearing the name "Bragg glass" for this property), which have 
indeed been observed in neutron diffraction experiments on Bi 2 Sr 2 CaCu20 8+ 5 (BSCCO) at 
low magnetic fields |TB[]. In the Bragg glass phase, the disordered FL array is modeled as an 



elastic manifold in a periodic random potential, similar to a randomly pinned charge density 
wave or a XY model in a random field [|l4| - [T6|Jl2]| . To give a thermodynamically stable phase, 
this requires the persistence of the topological order, or absence of unbound dislocations, even 
in the presence of disorder. In the neutron diffraction experiments by Cubitt et al. [13|], it has 



been observed that upon increasing the magnetic field, the Bragg peaks vanish, indicating 
an instability of the Bragg glass phase. Critical current measurements of Khaykovich et 
al. [17| show a sharp drop in the (local) critical current j c upon decreasing the magnetic 
induction below a critical value. This can be attributed to the "disentanglement" of FLs in 
the absence of dislocations when topological order is regained and dislocation loops vanish 
upon lowering the magnetic field. Similarly, changes in the I-V characteristics of YBCO 
based superlattices below a critical magnetic field can be interpreted as stemming from 
a sharp drop of the pinning energy and indicate a restoration of positional order in the 



VG |L8|. The existence of a topological transition has also been demonstrated in recent 



numerical studies fT9|,|20|j. In the closely related 3D XY model in a random field, vortex 



loops occur in a topological phase transition at a critical strength of the random field [19 



In simulations of disordered FL arrays |20| , a proliferation of dislocation lines has been found 



at a critical magnetic field in good agreement with the experimental results in Ref. [TB[ . 

Recently, the quantitative aspects of this issue have been addressed also analyti- 
cally [f2l| -|25|l- In Ref. |21] a self-consistent variational calculation and a scaling argument 



are presented, which show the topological stability of the elastic Bragg glass phase over a 
finite range of parameters that can be estimated by a Lindemann-like criterion (H). A more 
detailed discussion of the methods used in Ref. |2l] and their limitations as well as of the 



Lindemann-like criterion (|l|), which provides the basis for the calculations in the present 
work, will be given in the next paragraph. Recently, Fisher |22| has presented refined scal- 



ing arguments further supporting the existence of a topologically ordered Bragg glass phase. 
In Refs. p3|-p5[1, purely phenomenological Lindemann criteria are used as starting point 



for an estimate of the phase boundaries of the Bragg glass. Erta§ and Nelson [2~3| and 



Goldschmidt [^3J use "cage models" to mimic the interactions between FLs which yields an 
effective theory for a single FL in a random potential, to which they apply the conventional 
phenomenological Lindemann criterion. Giamarchi and Le Doussal apply a slightly 



modified phenomenological Lindemann criterion of the form (u 2 )(l) < c|Z 2 , where (u 2 )(l) is 
the (disorder-averaged) relative mean square displacement of neighbouring FLs separated by 
the FL spacing / (c^ is the Lindemann-number). However, it is well known that phenomeno- 
logical Lindemann criteria as used in Refs. |23H25| do not allow a theoretical description of 
a phase transition but can only give estimates of the location of the transition. They rely on 
the assumption that the phase transition reflects in the short scale behaviour of the system. 
Also the variational calculation and the scaling argument presented in Ref. |21| cannot give 
a complete description of the transition as only a detailed renormalization group (RG) anal- 
ysis of the problem would allow which is not yet available. The refined scaling arguments of 



Ref. 22 1 represent a further step towards this goal. 



In Ref. pT|] , the stability of the elastic Bragg glass phase has been investigated for a 
layered uniaxial geometry, where the magnetic field is parallel to the CuO-planes, by means 
of a self-consistent variational calculation and a scaling argument identifying the shear in- 
stability due to proliferating dislocations by the disorder-induced decoupling of the layers. 
For this geometry a Lindemann-like criterion has been derived, which is given below in ([!]) 
and relates the stability of the Bragg glass phase to the ratio of the positional correlation 
length and the FL spacing. These findings are supported by a more rigorous RG analysis 
|26| , |27H for a simplified model with only two layers of FLs in a parallel magnetic field. The 
usual experimental situation with the magnetic field perpendicular to the CuO-layers, which 
will be considered in the present work, is theoretically less understood mainly because dis- 
placements of the FLs have two components (biaxial) instead of one component (uniaxial). 
It is unclear at present whether topological phase transitions of the biaxial and uniaxial 
model share the same universality class p2 |. In Ref. [21J, it has been argued that the scaling 



argument for the uniaxial geometry can be generalized to the full, biaxial model leading to 
the following criterion estimating the range of stability of the Bragg glass phase with respect 
to a spontaneous formation of dislocation loops: 



I is the FL distance and A the magnetic penetration depth (we consider a magnetic field 
perpendicular to the CuO-planes of the HTSC, and will specify A below for such a geometry). 
Ri is the (transversal) positional correlation length of the disordered FL array, which is 
defined as the crossover length to the asymptotic large scale behaviour of the Bragg glass 
phase, where the average FL displacement starts to exceed the FL spacing /, see (^). c is 
a number, which was obtained in Ref. |^TJ to be of the order of c ~ (9(50), and ( « 1/5 is 



the roughness exponent of the pre-asymptotic so-called "random manifold" regime, see 
below. At the boundaries of the regime given by ([[]), a topological transition occurs, and 
dislocations proliferate. Beyond the transition line the FL array may form an amorphous 
VG with vanishing shear modulus or a viscous FL liquid. 



This article is divided into three parts. First, we will review the pre- asymptotic regimes 
of the FL array subject to point disorder on scales smaller than the positional correlation 
length Ri. This allows us to express Ri in terms of the microscopic parameters of the 
HTSC and the disorder strength, and to obtain its dependence on magnetic induction B 
and temperature T. In the second part, we will demonstrate the equivalence of the above 



criterion (jH) to the phenomenological Lindemann criterion in the form (u 2 )(l) < c 2 L l 2 (as for 



example used in Refs. p5|), where (u 2 )(l) is the (disorder-averaged) relative mean square 
displacement of neighboured FLs. This yields a relation c ~ c^ 2 between c from (|l|) and 
the Lindemann- number cl, and the value c ~ O(50) found by a variational calculation 



in Ref. [^TJ turns out to be in good agreement with a value cl ~ 0.15 widely used in 
the literature for the Lindemann-number. This equivalence further supports a scenario 
where the topological transition of the FL array subject to point disorder may be described 



as disorder- induced melting by unbound dislocations on the shortest scale / ||21|| . Finally, 
and most importantly from the experimental point of view, we estimate the region of the 
phase diagram of YBCO in the B-T plane (see Fig. 1) where the Bragg glass phase is 
stable and should be observable experimentally or numerically according to the Lindemann- 
like criterion ([!]). We find qualitative agreement with experiments The upper phase 
boundary of the Bragg glass, which we obtain using ([!]), turns out to be identical to the one 



obtained by Erta§ and Nelson [23 



II. POSITIONAL CORRELATION LENGTH 

To relate the positional correlation length Ri to the microscopic parameters of the 
HTSC and the disorder strength, we have to review the crossover between the different 
pre-asymptotic regimes of the dislocation-free disordered FL array preceding the asymp- 
totic Bragg glass phase, and the associated crossover length scales These crossovers 
are induced by the interplay between the FL interaction, the periodicity of the FL lattice 
and the disorder potential, which are in addition affected by thermal fluctuations, and lead 
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to essentially two different pre-asymptotic regimes: On the shortest scales, we have the 
"Larkin" or "random force" regime of Larkin and Ovchinnikov 0, which crosses over to 
the so-called "random manifold" regime at the Larkin length, before the asymptotic Bragg 
glass behaviour sets in on the largest scales exceeding the positional correlation length. In 
between this sequence of crossovers, one additional length scale is set by the FL interaction, 
which describes a crossover from a "single vortex" behaviour to a "collective" behaviour. 

In the following, we consider the usual experimental situation H||c of a magnetic field 
perpendicular to the CuO-planes of the HTSC FL positions are parameterized by the 
two-component displacement-field u(R, z) = u(r) in a continuum approximation of the 
Abrikosov lattice, where R is the vectors in the ab-plane and z is the coordinate in the 
c-direction, or by the Fourier transform u(K, k z ) = u(k). Let us adopt the convention 
to denote scales longitudinal to the FLs in the z-direction by L and transversal scales in 
the ab-plane by R. Moreover, it turns out to be convenient to use the reduced induction 
b = B/B c2 (T) = 2n^1 b /l 2 to measure the strength of the magnetic field. 



A. Interaction-induced Length Scale L* 

The dislocation- and disorder-free FL array can be described by elasticity theory (see 
Ref. I] for a review) in the displacement field u with the elastic moduli Cn, c 44 and c 66 , 
which can in general be dispersive (i.e., k-dependent in Fourier-space) due to the non- locality 
of the FL interaction. Except for extremely low magnetic fields, the FL lattice is essentially 
incompressible (cn ^> cqq), and we can neglect longitudinal compression modes to a good 
approximation. Note also, that the shear modulus Cqq is non-dispersive, because volume- 
preserving shear modes are not affected by the non-locality in the FL interaction. Then, the 
elastic Hamiltonian in the remaining transversal part ut of the displacement field is of the 
form: 

H ^t] = \l t^I^ {c 66 (K x u T ) 2 + c u [K, k z ](k z u T ) 2 } (2) 
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The dispersion- free shear modulus is given by cqq ~ eo/4/ 2 in the dense limit l/\ a b = 
(&/27t) _1 / 2 /k < 1 (with k = Aab/£ab), an d exponentially decaying c 66 oc exp (— //A a b)e /Z 2 in 
the dilute limit l/X a b > 1- eo = (^ > o/47rA a ( ) ) 2 is the basic energy (per length) scale of the FL. 
As estimates for YBCO we use £ o6 (0) « 15 A, e (0)^(0) » 1300K and k « 100 §. 

The tilt modulus c 44 = c 44 [i^, = c| 4 [i^] + c s u [k z ] is dispersive with the bulk- 
contribution 

C44[KJ -/ 2 1 + K 2 A 2 + A; 2 A 2 6 W 

dominating in the dense limit well within the Brillouin zone (BZ) K < Kbz — 
(approximated by a circular BZ) and the single vortex tilt modulus c 44 = c 44 [/c 2 ] [ p8| 

c 44 [^] = c 44 4" c 44 [k z ] 



-|^ 2 + xpf ln ( 1+A ^)J ( 4 ) 

dominating in the dilute limit and at the BZ boundaries K ~ Kbz or on scales R ~ /. 
£ — A a b/A c is the anisotropy ratio of the HTSC and approximately e ~ 1/5 in YBCO [|TJ. 
The single vortex tilt modulus has a strongly dispersive contribution c 44 em [/c 2 ] from the elec- 
tromagnetic coupling and an essentially dispersion-free contribution c 44 J from the Josephson 
coupling (where we neglect a logarithmically dispersive factor in c 44 J , which is of the order 
unity for the relevant wavevectors k z and magnetic inductions b). The length scale for the 
onset of dispersion in the bulk contribution c\ 4 is A c because elements of tilted FLs lying in 
the ab-plane will start to interact on scales R < A c [p9 ]. As length scale for the onset of 
dispersion, A c occurs as well in the Lindemann criterion ([I]). The contribution from the elec- 
tromagnetic coupling to the single vortex tilt modulus gives the local result c| 4 [0] ~ eo/l 2 for 
k z < l/Xab, but its strong dispersion c 44 m [/cj oc k~ 2 for k z > l/X a b leads to its suppression 
at small wavelengths k z > l/e\ a b where c 44 ~ c 44 J . 

From the competition of tilt and shear energy in (0), we can obtain a scaling relation 
between scales L longitudinal to the FLs and transversal scales R for typical fluctuations 
involving elastic deformation: 



^Jsjwar (5) 

V c 66 / 

The three-dimensional elastic Hamiltonian (|2]) is valid only on scales R > I or 

i>L ..f<i!V£l) 1/2 . (6) 

When we consider fluctuations on scales L < L* or R < I, the FL array breaks up into 
single FLs described by 1-dimensional elasticity in the longitudinal coordinate z with a line 
stiffness ei[k z ] = c s 44 [k z ]l 2 , because the shear energy containing the effects of FL interactions 
is always small compared to the tilt energy of the single FL. Thus the interaction-induced 
length scale L* separates a regime of "collective" behaviour described by 3D elasticity from 
a "single vortex" behaviour described by ID elasticity. L* starts to increase exponentially 
in the dilute limit l/X a b > 1 due to the exponential decay of c 66 . For the length scale 
L* given by ([]), we use therefore the local result c s 44 ~ c 44 m ~ eo/l 2 determined by the 
electromagnetic coupling. In the dense limit l/X a b < 1, the scale L* is smaller than e\ a b, 
and c s 44 ~ c 44 ~ e e 2 /l 2 , i.e., the dispersion-free contribution from the Josephson coupling 
dominates. This yields 

{/ < X ab : el 
l > : ( (V) exp 

for the interaction induced length scale L* in the dense and dilute limits. As we will show 
below, the criterion ([]]) is indeed equivalent to a Lindemann criterion in a more conventional 



form where fluctuations (u 2 ) on the transversal scale R ~ I are considered, see (133). There 



fore, the topological phase transition can be detected by considering fluctuations of single 
FLs on the longitudinal scale L ~ L*. 

We focus in this article on the upper branch of the topological transition line in moder- 
ately anisotropic compounds as YBCO such that the electromagnetic coupling can essentially 
be neglected. However, similar to the findings for thermal melting |3(J , the electromagnetic 
coupling and its strongly dispersive contribution to the single vortex tilt modulus plays an 
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important role for the disorder- induced topological phase transition in very anisotropic com- 
pounds such as BSCCO ||31||. We will consider effects from the electromagnetic coupling in 



detail in Ref. [p2 ]. We mention here only that in moderately anisotropic HTSC compounds 
with the upper branch bt, u {T) of the topological phase transition line will lie entirely in the 
dense regime b > 2ti/k 2 , but below the so-called "crossover field" b cr ~ {2ti / K 2 )(e\ 2 ab / d 2 ) 
above which L* < d and the layered structure of the HTSC becomes relevant at the tran- 
sition line and requires a discrete description in the c-direction. Only in this regime of 
magnetic inductions, the strongly dispersive electromagnetic contribution can be neglected 
at the topological transition (because L* < e\ a b), while a continuous description in the 
c-direction still applies. In the very anisotropic Bi-compounds, however, the upper branch 
b tiU (T) of the topological phase transition line typically lies in the dilute limit b < 2tt/k, 2 
where the electromagnetic coupling gives the relevant, strongly dispersive contribution to 
the single vortex tilt modulus c\ A ~ c^ m [k z ] oc k~ 2 . Because of this dispersion, the be- 
haviour of a single vortex of length L* changes drastically, and short-scale fluctuations on 
the (longitudinal) scale L ~ max{eX a b,d} give the main contribution [J32| . In the following, 
we consider the dense regime of a moderately anisotropic compounds such as YBCO and 
can thus neglect the dispersive electromagnetic contribution and use the dispersion-free, 
anisotropic result c s 44 ~ c^l ~ e^e 2 jl 2 . 

At the lower branch of the topological transition line in the dilute limit (where L* > X a b), 
we have to take into account the electromagnetic coupling and use the isotropic contribution 
c| 4 ~ C4 4 em [0] ~ eo 1 1 2 in the local limit. Furthermore, also in this regime effects from the 
strong dispersion of c^ m [k z ] have to be considered. The details of the calculation of the 



lower branch of the topological transition line will be given in Ref. ||32|| , and we will mention 
only the main results below. 
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B. Larkin Length 



When point disorder is introduced, every vortex at position H u in the Abrikosov lattice 
experiences a pinning potential V(r) with mean zero and short-range correlations 



V(r)V(r') = 7 *X.(R - - z>) , (8) 

where the overbar denotes an average over the quenched disorder. The strength of the 
disorder potential is given by 7 = n pin f pin1 where n pin is the density of pinning centers and 
f P i n the maximum pinning force exerted by one pinning center, and the effective range of 
the disorder potential is given by 

r T =(e ab + (u 2 ) th (0,h)y /2 (9) 

((. . .)th denotes a purely thermal average for a fixed ^(r)), which is equal to the size £ a & of 
the core of a vortex at T = but broadened by thermal fluctuations at higher temperatures. 
As proposed in Ref. [ffl], we introduce the dimensionless disorder strength 5 as 



it 



s = m- <io) 

The interaction with the disorder is described by the Hamiltonian 

HdisM = [ dzV(R„ + u(R„,z),z) . (11) 

For mean square displacements 

u(R, L) = ((u(r+(R,L))-u(r))2) 1/2 (12) 

smaller than the effective scale for variations of the disorder potential V, the FLs explore 
only one minimum of the disorder potential and perturbation theory in the displacements 
is valid. Expanding in ( |TT| ) the disorder potential V in u yields the random force theory of 
Larkin and Ovchinnikov ||. The (longitudinal) Larkin length is defined as the crossover 
scale for the random force regime, at which the average FL displacement becomes of order 
of the effective range r? of the point disorder: 
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(0,L e )~r T . (13) 



u 



It is important to note that for HTSCs such as YBCO and BSCCO, the generic disorder 
strength is such that 

L 5 < L* (14) 

in the range of magnetic inductions where the elastic Bragg glass will turn out to be sta- 
ble [^3j. Therefore, the random force regime lies entirely in the single vortex regime defined 
above, and the Larkin length is given by the single vortex result L|, which is at low 
temperatures [0] 

This result holds as long as r T ~ £ ab . However, above the depinning temperature TJ p of 
the single vortex, r T grows beyond ^ Q: 

4^d(l + exp((T/T^) 3 )) , (16) 

where the depinning temperature Tj is given by [0J 

TSp * ee U ^ ~ ee U 7 . (17) 

Above Tjp, L|(T) increases exponentially with temperature due to the fact that random 
forces are only marginally relevant for a single FL with two-component displacements 0: 

L\{T) * L|(0) \ T ^ ^ ' ) , , , (18) 

T > T| p : (T/TJ p ) exp ((T/Tj p ) J 

Let us discuss estimates of the quantities and T^ p at this point, which provide alternative 
measures of the disorder strength for a HTSC In Ref. ||23|| , the disorder strength is given by 
Ti p « 10K in BSCCO (where e w 1/100), which leads to 5/e « 1 with tfTfy. This estimate 
is considerably higher than typical values given in Ref. [|TJ for weak pinning. Therefore, we 
will use instead estimates in the range 8/e « 10 -3 . . . 10 _1 for YBCO in accordance with 
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Ref. which yield « 20 . . . 65K and values of the order of L|(0) « 30 . . . 6A for the 
(longitudinal) Larkin length in YBCO. 

For higher disorder strengths the T = Larkin length L|(0) can become smaller than 
the layer spacing d. In YBCO, where d ~ 12A, this happens for quite strong disorder 
S/e 2 ■ 10 -2 . Then, each pancake-like FL element of length d is pinned individually and 
we enter a strong pinning regime. This requires a description of pinning at the scale d as 
the smallest physical length scale in the longitudinal direction. In other words, we have 



to consider the disorder-induced relative displacement r\ = (u 2 ) of two pancake vortices in 
adjacent layers. This has been calculated in Ref. |34| at T = by means of an Imry-Ma 



argument (see also p5j) with the result 



rrf(0) -e ^ln(rfV^(0)) ln \2^eJ ' (19) 

where we introduced the mean-square disorder energy U 2 := jd^ b of a line-segment of 
length L ~ d. The result fll9|) is valid for r d (0) > £ ab , i.e., if the relative displacement 
exceeds the correlation length of the disorder potential, which is the case just for d > £|(0). 
The equation ([Tj^) has to be solved self-consistently, but in the following we will use the 
estimate obtained in the zeroth iteration 

S ^k^im) 3/2 ■ (20) 

The exponent 3/2 = 2^(1,0), see below (|22|), can be interpreted as the exponent charac- 
terizing the end-to-end displacement of a rigid rod that can tilt in a random potential. On 
scales L|(0) < L < d, each pancake can be treated as such a rigid rod of length L. 

Because the pinning is strong, each pancake remains individually pinned in the presence 
of thermal fluctuations until the thermal energy T is greater than the typical pinning energy 
Up of each pancake. Therefore, the result ( pD| ) remains to a good approximation valid in 
the whole temperature range T < U p : r^(T) ~ rd(0). This can be checked in a variational 
calculation along the lines of Ref. [35]. For YBCO with a disorder strength S/e ~ 2 • 10 _1 



we find U p ~ 70K. 
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Although the thermally increased Larkin length L|(T) becomes equal to the layer spacing 
d already at a temperature 

T L ,= d * n p (l - < U p , (21) 

[in YBCO, we find T^^d — 45K for 5/e ~ 2 • lO -1 ], the crossover from the strong pinning on 
the scale d to collective pinning on the scale of the Larkin length L S ^(T) can happen only at 
T ~ Up, where the strongly pinned individual pancakes can thermally depin. This result can 
be obtained from Ref. ||35|| , where it is shown that perturbation theory gives only a locally 
stable solution in a variational treatment of the pinning problem for two pancake vortices in 
adjacent layers in the temperature range T L ^ =d < T < U p whereas the result (|20|) represents 
the globally stable solution. 

C. Positional Correlation Length Ri 

On scales exceeding the Larkin length L|, the FLs start to explore many minima of the 
disorder potential V. However, as long as u(R, L) is smaller than the FL spacing I, FLs are 
not competing for the same minima, and the FLs experience effectively independent disorder 
potentials. This leads to the approximation TC^is [ u ] ~ / d 3 rV(r, u(r)) (on longitudinal scales 
exceeding the layer spacing d), where V has also short-range correlations in u. This regime 
is referred to as the random manifold regime |12|]. For a d-dimensional (dispersion- free) 
elastic manifold with a n-component displacement field u, the scaling behaviour of the (uu)- 
correlations is known to be 

u(0, L) ~ L^ d ' n) (22) 

with a roughness exponent ((d, n). We are interested here in the case d — 1, n — 2, which is 
realized on scales d, Li < L < L* in the single vortex regime, where the FLs are described 
as 1-dimensional elastic manifolds, and the case d = 3, n = 2 on scales L* < L < Li (or 
transversal scales I < R < R{) in the collective regime, where the FL array is described as 
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3-dimensional elastic manifold. Li and Ri are the positional correlation lengths, which are 
defined as the crossover scales for the random manifold regime, at which the average FL 
displacement becomes of the order of the FL distance I: 

u(R l ,0) = u(0,L l ) = l . (23) 

On scales R > Ri, where u(R) > I, FLs start to compete for the same minima, and the 
periodicity of the FL lattice becomes crucial |5|Jl2]|. The FL array reaches its asymptotic 



behaviour of the Bragg glass phase with only logarithmically diverging (tm)-correlations, i.e., 
quasi long range positional order. The best estimates available for the roughness exponents 



are 36 



C(l,2)«5/8 and ( = C(3, 2) « 1/5 , (24) 

where the latter occurs also in the above Lindemann criterion ([!]). In the collective regime 
the scaling relation ( p2|) gets slightly modified by the dispersion (|3|) of c| 4 to 

>C(3,2) 



u 2 



(R,0)~(\l + R^> Z> , (25) 



as can be checked by means of a simple Flory-type argument, where we equate the typical 
disorder energy and elastic energy @ on one dominant scale. (As suggested by a more elab- 
orate variational calculation as in Ref. |12| we neglect possible small logarithmic corrections 
of order In (1/e) in fl25]).) Note that for I < R < A c , the relative displacements (|25|) are only 
marginally growing due to the dispersion of c\ A . 

The scaling relations (|22"| , |2"5"D enable us to obtain the relation between the (transversal) 



positional correlation length Ri and the (longitudinal) Larkin length L|, which will allow us 
to express Ri in terms of microscopic parameters, both for weak pinning on the scale L| (for 
L^(T) > d) and for strong pinning of pancakes on the scale d (for L S ^(T) < d). 



Applying the scaling relation ( |2"2"D for the (im)-correlations to the single vortex random 
manifold regime on longitudinal scales L| < L < L*, we obtain for the case £|(0) > d of 
weak pinning 

14 



/ \C(1,2) 

u* = 0) ~ u(0, L*) ~ r r (^^i^yj • (26) 

In the same manner we can use ( P5| ) in the collective random manifold regime on transversal 
scales I < R < Rf. 

( \2 i p2\ C(3,2) R 2C(3,2) 

f =" ! «)-'(w) ^-(f) (27) 

with Ri~^> \ c ~^> I. Using (|6|,^7]), Ri can be expressed as 

/ , v i/C(3,2) / rsmN c(i,2)/C(3,2) 

*<^(£) (^) • < 28 > 

With the results © for r T , © for L*, and ([18]) for L|(T) together with £(3,2) « 1/5 and 
C(l,2) ~ 5/8 fl24|), this yields the desired expression for i^: 

/ , \ -15/16 , x \ -25/24 

fr<T? : 1 

[ T > 7- : (T/Tj p ) exp (j (T/Tj p ) J 

The weakening of the pinning by thermal fluctuations leads to an exponential increase of 
Ri(T) for temperatures above the depinning temperature TL similar (and related) to the be- 
haviour of the thermally increased Larkin length L|(T). For inductions b/2ii = 10~ 4 . . . 10~ 2 
in the dense limit 6/27T > 1/k 2 , a disorder strength 5/e ~ 10~ 2 , and A c (0) ~ 7500A, 
we obtain (transversal) positional correlation lengths Ri(0) ~ (10 4 ...10 6 ) • A c ~ 7, 5 • 
(1CT 1 . . . 10)cm, which are extremely large indicating that over a wide range of length scales 
the pre-asymptotic random manifold regimes should be observable rather than the asymp- 
totic Bragg glass regime. 

For the case L|(0) < d of strong pinning of pancakes on the scale d, we apply ( p^ ) for 



the ('U'u)-correlations to the single vortex random manifold regime on longitudinal scales 
d < L < L* and obtain for low temperatures 

u* = u{R = I, 0) ~ u(0, L = L*) ~ r d (0) — (30) 



d 
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instead of (|26|). Using this and (|27D, Ri can be expressed as 



/ , v 1/C(3,2) / , x C(l,2)/C(3,2) 

*<°>^U W • (31) 



With the results fl20|) for r d (0), © for L*, and C(3,2) « 1/5 and £(1,2) w 5/8 we 
obtain for the positional correlation length Ri. 

f h \ ~ 15/16 /X\ ~ 5/6 ft \ 5/8 

This result is valid for temperatures T < U p and gives a temperature independent, smaller 
value than ( p9|) in this temperature range. At T ~ C/ p , pancakes can thermally depin for 
strong pinning, and we expect a crossover to the weak pinning result (^) with a pronounced 
increase of Ri{T) with temperature. In YBCO, strong pinning is realized for S/e > 2 • 10~ 2 . 
For 8/e w 2 • 10 _1 and with the layer spacing d w 12A, we find Ri(0) ~ 5 • (10 2 . . . 10 4 ) • A c ps 
3, 8 • (10 -2 . . . l)cm in the induction range b/2-K = 10 -4 . . . 10~ 2 in the dense limit. 

III. LINDEMANN CRITERION 

Let us now show the equivalence of the Lindemann-like criterion ([!]) obtained in Ref. [21 



to the conventional form of the Lindemann criterion generalized to a disordered system. The 
Lindemann criterion has been proven as a very efficient phenomenological tool to obtain the 
thermal melting curves of lattices, e.g. the disorder-free FL lattice. There, it is formulated 
in its conventional form 

(u 2 ) th = c\l 2 , (33) 

with a Lindemann-number cl ~ 0.1 . . .0.2. For the thermal melting of the FL array, the 
main contributions to the left hand side of ([!3|) come from fluctuations on the shortest scale, 
which is in the transverse direction the FL spacing I, i.e., (u 2 ) t h ~ (u 2 )th(l, 0) (note that we 
apply again a convention like fll2|)). Therefore, the straightforward generalization of ([33]) to 
the disorder- induced melting by dislocations is 
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( M 2)(/,0)^( M 2)(o,L*)=^ = ci/ 2 , (34) 



where we consider again fluctuations on the shortest scale R ~ I. In Ref. [pi]] , one derivation 



of the criterion (|1]) was based on a variational calculation for a layered superconductor in 
a parallel field. There it was found, that unbound dislocations proliferate indeed on the 
shortest scale at the topological transition described by ([!]), i.e., in between every layer and 
thus with a distance I. This suggests that the use of the conventional phenomenological 
Lindemann criterion in the form (]34] ) might be one possibility to obtain the topological 
transition line. 

This can be further justified by showing that the criterion ([!]), obtained in Ref. |^T] 



on 



the basis of a scaling argument and a variational calculation for a uniaxial model, is actually 



equivalent to the phenomenological Lindemann criterion fl3~4|): Considering the relation 
between u* and /, it becomes clear that (p]) is the analog of the Lindemann criterion (|34"D 
formulated in terms of the underlying transversal scales rather than the corresponding dis- 
placements. Using (|27|), the criterion (|I]) for the stability of the Bragg glass can be written 

as 

ul < c-H 2 . (35) 

This is just the above phenomenological Lindemann criterion (|3~4"1), and we can identify 

c^cl 2 . (36) 

We see that the equivalence of the criterion ([!]) to the phenomenological Lindemann cri- 
terion (|34]) includes the agreement of the appearing numerical factors: The value for the 
Lindemann-number cl ~ 0.15, widely used in the literature, produces a good agreement in 
(|36|) with the value c ~ O(50) obtained by the variational calculation. This equivalence to 
a scenario where disorder- induced fluctuations on the shortest scale "melt" the FL array 
favours a first order transition scenario for the topological transition, which could not be 
excluded in the experiments |T7j]. As we will see, the quantity ul is equivalent to the mean 
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square displacement of the "effective" FL studied in the "cage model" of Erta§ and Nel- 
son ||23|| . They apply the Lindemann criterion directly in its phenomenological form (|54"D to 
the "caged" FL. Using ( |31)| , |2"4"|) , we can cast the Lindemann-like criterion ([!]) into the form 

Rt > cl lK (l 2 + A c 2 ) 1/2 « max {/, AJ . (37) 



IV. PHASE DIAGRAM 

Let us now address the issue of phase boundaries of the topologically ordered Bragg glass 
in the B-T plane as they follow from the Lindemann-like criterion ([J) in the above form ([37|) . 
The results are summarized in Fig. 1. The boundary of the regime given by ( j37|) defines 
a topological transition line B t (T), where dislocations proliferate and the topological order 
of the Bragg glass phase is lost. The upper branch b ttU (T) of this line can be obtained by 
applying the expressions (E^) or fl32|) for the positional correlation length Ri in the dense 



limit b > 2tt/k to the criterion (|37|). 

For weak pinning or L|(0) > d such that we have collective pinning on the scale L|, this 
yields a condition b < b tiU (T) in the b-T plane with 



(?) 


-10/9 

16/3 
C L 




' T « 7£ 
J>TS P 


(0)< 



2tt 



ee ojab \ ^ 16/3 



, 10/3 exp (I (t/T^ 3 



(38) 



Note that the transition line (^) is identical to the one obtained by Erta§ and Nelson 
by applying the conventional phenomenological Lindemann criterion to a "cage model" for a 
single FL (this demonstrates the equivalence of the displacement as defined in ([33]) to the 
average displacement of the "caged" FL). Estimates of b ttU (0) strongly depend on the chosen 
value for the Lindemann number cl ~ 0.1 . . . 0.2. A value cl ~ 0.15 and a disorder strength 
5/e w 10^ 2 lead to 6t,«(0) w 4 ■ 10^ 2 or ^(0) » 6T with B c2 {0) w 150T. For temperatures 



T < TJp the transition line is essentially temperature-independent because the mechanism 



for the proliferation of dislocation loops is purely disorder-driven at low temperatures jzT 
For T > TJp it increases exponentially due to the very effective weakening of the pinning 
effects by thermal fluctuations in the single vortex regime. 

For L|(0) < d, i.e., strong pinning of pancakes on the scale d (realized for 8/e ^ 2 • 1CT 2 
in YBCO), we obtain instead 

K ^)^2^(^ il3 ^Y (39) 

at low temperatures. This result gives a lower induction for the topological transition than 
fl38|) and remains valid up to the temperature T ~ U p , where pancakes can thermally depin. 
For T > Up we expect a pronounced increase of b t , u (T) with temperature and a crossover to 
the weak pinning result Q38|), cf. Fig. 1. Estimates for b ttU (0) are again very susceptible to 
changes in the chosen value for the Lindemann number cl- For cl ~ 0.15 and 8/e ~ 2- 1CT 1 , 
we find b t>u {0) « 8 • 10~ 4 or 5^(0) « 0, 12T for YBCO. 

The estimates for B tjU (0) obtained from (|38|) and (|39|) are in qualitative agreement with 
the experimental results for the magnetic field where in YBCO based superlattices a change 
in the I-V characteristics has been observed |I8 |. Both (|38|) and (|39|) show that the magnetic 



induction bt jU at the topological transition decreases for stronger anisotropy or effectively 
larger disorder strength 8/e, and the stability region of the topologically ordered Bragg glass 
shrinks. 

At some temperature T XiU above T^ p (and T L(=d ), the topological transition line b t u {T) 
will terminate in the upper branch of the melting curve b mu (T), which is 

b m , u (T) « 2*4 (40) 



in this regime of inductions for the moderately anisotropic compound YBCO |30|. The 
temperature T x<tl can be determined from the condition that the thermally increased Larkin 
length L S ^(T) becomes equal to the scale L* of the dominant fluctuations at the melting 



line and the topological transition line. Because (u 2 )(0, L|(T)) = (■u 2 ) t / 1 (0, L|(T)) at the 
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thermally increased Larkin length, the Lindemann criteria ([33|) for thermal melting and 
(|37|) in the form ( |34| ) for the topological phase transition are indeed fulfilled simultaneously 
if L* = L|(T): 



(u*)(0,L*) = (u 2 ) th (0,L*) = c 2 L l 2 . (41) 



This yields 



'^ 2/3 



T^TMn 1 / 3 H c L 2 (42) 



In YBCO, we find T XjU « 80K with the above estimates c L « 0.15 and d>/e ~ 2 ■ 10 _1 . 
For T > T XjU beyond the melting curve b mM (T), the FL array melts into a disordered FL 
liquid, and the Bragg glass order is destroyed by the thermal fluctuations on small scales 
where disorder- induced fluctuations are irrelevant, whereas above the transition line b ttU (T) 
the Bragg glass "melts" by disorder-induced fluctuations, when unbound dislocation loops 
proliferate. For T < T x u , we find b m>u {T) > b tiU (T), and the melting curve b m ^ u {T) lies above 
the topological transition line in the b-T plane. Therefore, we expect for temperatures 
T < T X)U that an amorphous, i.e., topologically disordered vortex glass melts into a vortex 
liquid at the thermal melting line b myU (T) and consequently, that the melting transition into 
a vortex liquid at b m , u (T) is of a different nature below and above T x>u . In the experiments 
reported in Ref. ||11|| , such a change in the properties of the melting transition has indeed 
been observed in YBCO at a temperature around 75K which is in fairly good agreement 
with our result for T x ^ u . 

Let us now give the main results for the lower branch of the topological transition line 
bt } i(T) at which the strongly dispersive contribution from the electromagnetic coupling to the 
single vortex tilt modulus is dominating. At low inductions in the dilute limit b <C 2tt/k 2 , 
the criterion (|37| ) will be violated due to the exponential decrease of the shear modulus 
C66, or increase of the interaction-induced length scale L* (0). At low temperatures T ~ 0, 
the positional correlation length Ri(0) can be determined also from (EH) using the isotropic 



single vortex Larkin length (given by (|15|) with e = 1) and the appropriate result for L*. 
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Because the isotropic Larkin length is always greater than the layer spacing, the layered 
structure is irrelevant for the collective pinning at low inductions. The criterion ( |3TD then 
yields for the lower branch of the topological transition line 

6 M (0)«^ln- 2 (cf 5 K 6 / 5 5(0)- 2 / 3 ) (43) 

Due to the strong dispersion of c s ^ m [k z ], the thermal depinning at higher temperatures 
is more complex and involves several crossover temperatures. However, only above the 
isotropic single vortex depinning temperature TJ , (given by (|17D with e = 1) the positional 
correlation length is increasing exponentially similarly to the thermally increased isotropic 
Larkin length. This gives only a weak logarithmic temperature dependence for T < T| pi 
whereas we find the asymptotics 

(44, 

at temperatures T ^> T$ pi well above the isotropic single vortex depinning temperature. At 
a temperature T x j (> Td Pji ), b t j(T) will terminate in the lower branch of the melting curve 



b m ,i(T), which increases logarithmically with temperature pO] 



b mil (T) « ^ln- 2 (^|^) . (45) 



K 2 



Analogously to the findings for the upper branch of the topological transition line, T x j can be 
determined from the condition that the thermally increased isotropic Larkin length becomes 
equal to the scale L* at the melting line. This yields 

TLAJI* [ cW* ( 5 -\ I (JO) 



x,l — 1 dp,i 



£ 



With cl ~ 0.15 and 8/e w 10 2 , we obtain b t ,i(0) ~ 0.16(27t/k 2 ) « 1 ■ 10 4 , which is by a 
factor of 40 smaller than &t lU (0) and experimentally hard to verify due to the small inductions 
B t<u (0) w 150G. Furthermore we find T dP)i w 70K and T Xt i w 85K. From (|43|) it is clear that 
the transition line bt >u {T) increases with the disorder strength so that the stability region of 
the topologically ordered Bragg glass shrinks. 
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V. CONCLUSION 



In conclusion, we have obtained the region in the phase diagram of YBCO in the B-T 
plane, where the topologically ordered vortex glass should be observable, and the topological 
transition lines B tjU (T) and B t j(T), where dislocation loops proliferate. The resulting phase 
diagram, as given by the formulae (^), fl3D|), fl4"3]), and ([PJ) is depicted in Fig. 1. The results 



are in qualitative agreement with the experimental data of Ref. [ | if the observed changes 



in the I-V characteristics are attributed to a topological transition of the disordered vortex 
array. 

The phase diagram is based on the Lindemann-like criterion (P or (|37|) , which has been 
obtained by a variational calculation for a uniaxial model and a scaling argument presented 
in Ref. plR . We have demonstrated the equivalence to the conventional phenomenological 
formulation of the Lindemann criterion ([34]) up to the involved numerical factors, i.e., the 
Lindemann- number c^. Our results for the upper branch of the topological transition line 



B t ^ u {T) agree with Ref. |23|, where the conventional phenomenological Lindemann-criterion 



was applied to the disorder- induced "melting" in the framework of a "cage model" . 
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Fig. 1: Schematic phase diagram in the b-T plane {b = B/B c2 (T)) showing the stability 
regime of the topologically ordered Bragg glass phase (hatched region). Its phase boundaries 
are given by the upper and lower branch b ttU (T) and b t j(T) (solid lines) of a topological 
transition line, where dislocations proliferate. They terminate in the two branches b mtU (T) 
and b m j(T) (dashed lines) of the melting curve, where the FL array melts by thermal 
fluctuations into a (disordered) FL liquid. 



26 



